A mutation associated with resistance to synthetic pyrethroids is widespread in US populations of the tropical lineage of Rhipicephalus sanguineus s.l.

The brown dog tick, Rhipicephalus sanguineus sensu lato (s.l.), is an important vector for Rickettsia rickettsii, causative agent of Rocky Mountain spotted fever. Current public health prevention and control efforts to protect people involve preventing tick infestations on domestic animals and in and around houses. Primary prevention tools rely on acaricides, often synthetic pyrethroids (SPs); resistance to this chemical class is widespread in ticks and other arthropods. Rhipicephalus sanguineus s.l. is a complex that likely contains multiple unique species and although the distribution of this complex is global, there are differences in morphology, ecology, and perhaps vector competence among these major lineages. Two major lineages within Rh. sanguineus s.l., commonly referred to as temperate and tropical, have been documented from multiple locations in North America, but are thought to occupy different ecological niches. To evaluate potential acaricide resistance and better define the distributions of the tropical and temperate lineages throughout the US and in northern Mexico, we employed a highly multiplexed amplicon sequencing approach to characterize sequence diversity at: 1) three loci within the voltage-gated sodium channel (VGSC) gene, which contains numerous genetic mutations associated with resistance to SPs; 2) a region of the gamma-aminobutyric acid-gated chloride channel gene (GABA-Cl) containing several mutations associated with dieldrin/fipronil resistance in other species; and 3) three mitochondrial genes (COI, 12S, and 16S). We utilized a geographically diverse set of Rh sanguineus s.l. collected from domestic pets in the US in 2013 and a smaller set of ticks collected from canines in Baja California, Mexico in 2021. We determined that a single nucleotide polymorphism (T2134C) in domain III segment 6 of the VGSC, which has previously been associated with SP resistance in Rh. sanguineus s.l., was widespread and abundant in tropical lineage ticks (>50 %) but absent from the temperate lineage, suggesting that resistance to SPs may be common in the tropical lineage. We found evidence of multiple copies of GABA-Cl in ticks from both lineages, with some copies containing mutations associated with fipronil resistance in other species, but the effects of these patterns on fipronil resistance in Rh. sanguineus s.l. are currently unknown. The tropical lineage was abundant and geographically widespread, accounting for 79 % of analyzed ticks and present at 13/14 collection sites. The temperate and tropical lineages co-occurred in four US states, and as far north as New York. None of the ticks we examined were positive for Rickettsia rickettsii or Rickettsia massiliae.

Pan-Enterovirus Characterization Reveals Cryptic Circulation of Clinically Relevant Subtypes in Arizona Wastewater.

Background: Most seasonally circulating enteroviruses result in asymptomatic or mildly symptomatic infections. In rare cases, however, infection with some subtypes can result in paralysis or death. Of the 300 subtypes known, only poliovirus is reportable, limiting our understanding of the distribution of other enteroviruses that can cause clinical disease. Objective: The overarching objectives of this study were to: 1) describe the distribution of enteroviruses in Arizona during the late summer and fall of 2022, the time of year when they are thought to be most abundant, and 2) demonstrate the utility of viral pan-assay approaches for semi-agnostic discovery that can be followed up by more targeted assays and phylogenomics. Methods: This study utilizes pooled nasal samples collected from school-aged children and long-term care facility residents, and wastewater from multiple locations in Arizona during July-October of 2022. We used PCR to amplify and sequence a region common to all enteroviruses, followed by species-level bioinformatic characterization using the QIIME 2 platform. For Enterovirus-D68 (EV-D68), detection was carried out using RT-qPCR, followed by confirmation using near-complete whole EV-D68 genome sequencing using a newly designed tiled amplicon approach. Results: In the late summer and early fall of 2022, multiple enterovirus species were identified in Arizona wastewater, with Coxsackievirus A6, EV-D68, and Coxsackievirus A19 composing 86% of the characterized reads sequenced. While EV-D68 was not identified in pooled human nasal samples, and the only reported acute flaccid myelitis case in Arizona did not test positive for the virus, an in-depth analysis of EV-D68 in wastewater revealed that the virus was circulating from August through mid-October. A phylogenetic analysis on this relatively limited dataset revealed just a few importations into the state, with a single clade indicating local circulation. Significance: This study further supports the utility of wastewater-based epidemiology to identify potential public health threats. Our further investigations into EV-D68 shows how these data might help inform healthcare diagnoses for children presenting with concerning neurological symptoms.

Large Clusters of Invasive emm49 Group A Streptococcus Identified within and across Arizona Healthcare Facilities through Statewide Genomic Surveillance System, 2019-2021.

A statewide genomic surveillance system for invasive Group A Streptococcus was implemented in Arizona in June 2019, resulting in 1,046 isolates being submitted for genomic analysis to characterize emm-types and identify transmission clusters. Eleven of the 32 identified distinct emm-types comprised >80% of samples, with 29.7% of all isolates being typed as emm49 (and its genetic derivative emm151). Phylogenetic analysis initially identified an emm49 genomic cluster of four isolates that rapidly expanded over subsequent months (June 2019-February 2020). Public health investigations identified epidemiologic links with three different long-term care facilities, resulting in specific interventions. Unbiased genomic surveillance allowed for identification and response to clusters that would have otherwise remained undetected.

Bedaquiline resistance in patients with drug-resistant tuberculosis in Cape Town, South Africa: a retrospective longitudinal cohort study.

BACKGROUND: Bedaquiline is a life-saving tuberculosis drug undergoing global scale-up. People at risk of weak tuberculosis drug regimens are a priority for novel drug access despite the potential source of Mycobacterium tuberculosis-resistant strains. We aimed to characterise bedaquiline resistance in individuals who had sustained culture positivity during bedaquiline-based treatment. METHODS: We did a retrospective longitudinal cohort study of adults (aged ≥18 years) with culture-positive pulmonary tuberculosis who received at least 4 months of a bedaquiline-containing regimen from 12 drug-resistant tuberculosis treatment facilities in Cape Town, South Africa, between Jan 20, 2016, and Nov 20, 2017. Sputum was programmatically collected at baseline (ie, before bedaquiline initiation) and each month to monitor treatment response per the national algorithm. The last available isolate from the sputum collected at or after 4 months of bedaquiline was designated the follow-up isolate. Phenotypic drug susceptibility testing for bedaquiline was done on baseline and follow-up isolates in MGIT960 media (WHO-recommended critical concentration of 1 µg/mL). Targeted deep sequencing for Rv0678, atpE, and pepQ, as well as whole-genome sequencing were also done. FINDINGS: In total, 40 (31%) of 129 patients from an estimated pool were eligible for this study. Overall, three (8%) of 38 patients assessable by phenotypic drug susceptibility testing for bedaquiline had primary resistance, 18 (47%) gained resistance (acquired or reinfection), and 17 (45%) were susceptible at both baseline and follow-up. Several Rv0678 and pepQ single-nucleotide polymorphisms and indels were associated with resistance. Although variants occurred in Rv0676c and Rv1979c, these variants were not associated with resistance. Targeted deep sequencing detected low-level variants undetected by whole-genome sequencing; however, none were in genes without variants already detected by whole-genome sequencing. Patients with baseline fluoroquinolone resistance, clofazimine exposure, and four or less effective drugs were more likely to have bedaquiline-resistant gain. Resistance gain was primarily due to acquisition; however, some reinfection by resistant strains occurred. INTERPRETATION: Bedaquiline-resistance gain, for which we identified risk factors, was common in these programmatically treated patients with sustained culture positivity. Our study highlights risks associated with implementing life-saving new drugs and shows evidence of bedaquiline-resistance transmission. Routine drug susceptibility testing should urgently accompany scale-up of new drugs; however, rapid drug susceptibility testing for bedaquiline remains challenging given the diversity of variants observed. FUNDING: Doris Duke Charitable Foundation, US National Institute of Allergy and Infectious Diseases, South African Medical Research Council, National Research Foundation, Research Foundation Flanders, Stellenbosch University Faculty of Medicine Health Sciences, South African National Research Foundation, Swiss National Science Foundation, and Wellcome Trust.

Unique Genomic Epidemiology of COVID-19 in the White Mountain Apache Tribe, April to August 2020, Arizona.

The first case of coronavirus disease 2019 (COVID-19) within the White Mountain Apache Tribe (WMAT) in Arizona was diagnosed almost 1 month after community transmission was recognized in the state. Aggressive contact tracing allowed for robust genomic epidemiology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and subsequent phylogenetic analyses implicated only two virus introductions, which resulted in the spread of two unique viral lineages on the reservation. The phylogenies of these lineages reflect the nature of the introductions, the remoteness of the community, and the extraordinarily high attack rates. The timing and space-limited nature of the outbreaks validate the public health tracing efforts involved, which were illustrated by multiple short transmission chains over a period of several weeks, eventually resulting in extinction of the lineages. Comprehensive sampling and successful infection control efforts are illustrated in both the effective population size analyses and the limited mortality outcomes. The rapid spread and high attack rates of the two lineages may be due to a combination of sociological determinants of the WMAT and a seemingly enhanced transmissibility. The SARS-CoV-2 genomic epidemiology of the WMAT demonstrates a unique local history of the pandemic and highlights the extraordinary and successful efforts of their public health response. IMPORTANCE This article discusses the introduction and spread of two unique viral lineages of SARS-CoV-2 within the White Mountain Apache Tribe in Arizona. Both genomic sequencing and traditional epidemiological strategies (e.g., contract tracing) were used to understand the nature of the spread of both lineages. Beyond providing a robust genomic analysis of the epidemiology of the outbreaks, this work also highlights the successful efforts of the local public health response.

St. Louis Encephalitis Virus in the Southwestern United States: A Phylogeographic Case for a Multi-Variant Introduction Event.

Since the reemergence of St. Louis Encephalitis (SLE) Virus (SLEV) in the Southwest United States, identified during the 2015 outbreak in Arizona, SLEV has been seasonally detected within Culex spp. populations throughout the Southwest United States. Previous work revealed the 2015 outbreak was caused by an importation of SLEV genotype III, which had only been detected previously in Argentina. However, little is known about when the importation occurred or the transmission and genetic dynamics since its arrival into the Southwest. In this study, we sought to determine whether the annual detection of SLEV in the Southwest is due to enzootic cycling or new importations. To address this question, we analyzed 174 SLEV genomes (142 sequenced as part of this study) using Bayesian phylogenetic analyses to estimate the date of arrival into the American Southwest and characterize the underlying population structure of SLEV. Phylogenetic clustering showed that SLEV variants circulating in Maricopa and Riverside counties form two distinct populations with little evidence of inter-county transmission since the onset of the outbreak. Alternatively, it appears that in 2019, Yuma and Clark counties experienced annual importations of SLEV that originated in Riverside and Maricopa counties. Finally, the earliest representatives of SLEV genotype III in the Southwest form a polytomy that includes both California and Arizona samples. We propose that the initial outbreak most likely resulted from the importation of a population of SLEV genotype III variants, perhaps in multiple birds, possibly multiple species, migrating north in 2013, rather than a single variant introduced by one bird.

Applying Genomic Epidemiology to Characterize a COVID-19 Outbreak in a Developmentally Disabled Adult Group Home Setting, Arizona.

Individuals living in congregate settings, including those in group homes, have been disproportionately impacted by COVID-19 and may be at increased risk of exposure or infection due to underlying illness. In mid-May 2020, local public health officials responded to an outbreak of COVID-19 among staff and residents associated with a multi-residential group home that provides care for adults with intellectual and developmental disabilities. Samples were collected at 16 of the homes. In four of the homes all the residents tested positive, and in the remaining 12 houses where samples were collected, all residents tested negative. Of the 152 individuals tested, 15/58 (25.9%) residents and 27/94 (28.7%) staff were positive for SARS-CoV-2, including eight hospitalizations and four deaths. Phylogenetic analysis of genomes from this outbreak in the context of genomes from Northern Arizona shows that very few mutations separate the samples from this outbreak. A potential transmission network was developed to illustrate person-place epidemiologic linkages and further demonstrates the dynamic connections between staff and residents with respect to each group home location. Epidemiologic and genomic evidence correlate, and suggest that asymptomatic infected staff likely introduced and spread COVID-19 in this setting. Implementation of public health prevention measures alongside rapid genomic analysis can help guide policy development and guide management efforts to prevent and mitigate future outbreaks.

Methods for sequencing the pandemic: benefits of rapid or high-throughput processing.

Genomic epidemiology has proven successful for real-time and retrospective monitoring of small and large-scale outbreaks. Here, we report two genomic sequencing and analysis strategies for rapid-turnaround or high-throughput processing of metagenomic samples. The rapid-turnaround method was designed to provide a quick phylogenetic snapshot of samples at the heart of active outbreaks, and has a total turnaround time of <48 hours from raw sample to analyzed data. The high-throughput method was designed for semi-retrospective data analysis, and is both cost effective and highly scalable. Though these methods were developed and utilized for the SARS-CoV-2 pandemic response in Arizona, U.S, and we envision their use for infectious disease epidemiology in the 21 st Century.

Novel Odoribacter splanchnicus Strain and Its Outer Membrane Vesicles Exert Immunoregulatory Effects in vitro.

Odoribacter splanchnicus, belonging to the order Bacteroidales, is a common, short-chain fatty acid producing member of the human intestinal microbiota. A decreased abundance of Odoribacter has been linked to different microbiota-associated diseases, such as non-alcoholic fatty liver disease, cystic fibrosis and inflammatory bowel disease (IBD). The type strain of O. splanchnicus has been genome-sequenced, but otherwise very little is known about this anaerobic bacterium. The species surfaces in many microbiota studies and, consequently, comprehension on its interactions with the host is needed. In this study, we isolated a novel strain of O. splanchnicus from a healthy fecal donor, identified it by genome sequencing and addressed its adhesive, epithelium reinforcing and immunoregulatory properties. Our results show that O. splanchnicus strain 57 is non-adherent to enterocytes or mucus, does not reinforce nor compromise Caco-2 monolayer integrity and most likely harbors penta-acylated, less endotoxic lipid A as part of its lipopolysaccharide (LPS) structure based on the lack of gene lpxM and in vitro results on low-level NF-κB activity. The studies by transmission electron microscopy revealed that O. splanchnicus produces outer membrane vesicles (OMV). O. splanchnicus cells, culture supernatant i.e., spent medium or OMVs did not induce interleukin-8 (IL-8) response in HT-29 enterocyte cells suggesting a very low proinflammatory capacity. On the contrary, the treatment of HT-29 cells with O. splanchnicus cells, spent medium or OMVs prior to exposure to Escherichia coli LPS elicited a significant decrease in IL-8 production as compared to E. coli LPS treatment alone. Moreover, O. splanchnicus spent supernatant induced IL-10 production by immune cells, suggesting anti-inflammatory activity. Our in vitro findings indicate that O. splanchnicus and its effector molecules transported in OMVs could potentially exert anti-inflammatory action in the gut epithelium. Taken together, O. splanchnicus seems to be a commensal with a primarily beneficial interaction with the host.

Molecular Evaluation of Fluoroquinolone Resistance in Serial Mycobacterium tuberculosis Isolates from Individuals Diagnosed with Multidrug-Resistant Tuberculosis.

Fluoroquinolones (FQ) are crucial components of multidrug-resistant tuberculosis (MDR TB) treatment. Differing levels of resistance are associated with specific mutations within the quinolone-resistance-determining region (QRDR) of gyrA We sequenced the QRDR from serial isolates of MDR TB patients in the Preserving Effective TB Treatment Study (PETTS) with baseline FQ resistance (FQR) or acquired FQ resistance (FQACQR) using an Ion Torrent Personal Genome Machine (PGM) to a depth of 10,000× and reported single nucleotide polymorphisms in ≥1% of reads. FQR isolates harbored 15 distinct alleles with 1.3 (maximum = 6) on average per isolate. Eighteen alleles were identified in FQACQR isolates with an average of 1.6 (maximum = 9) per isolate. Isolates from 78% of FQACQR individuals had mutant alleles identified within 6 months of treatment initiation. Asp94Gly was the predominant allele in the initial FQ-resistant isolates followed by Ala90Val. Seventy-seven percent (36/47) of FQACQR group patients had isolates with FQ resistance alleles prior to changes to the FQ component of their treatment. Unlike the individuals treated initially with other FQs, none of the 21 individuals treated initially with levofloxacin developed genotypic or phenotypic FQ resistance, although country of residence was likely a contributing factor since 69% of these individuals were from a single country. Initial detection of phenotypic resistance and genotypic resistance occurred simultaneously for most; however, phenotypic resistance occurred earlier in isolates harboring mixtures of alleles of very low abundance (<1% of reads), whereas genotypic resistance often occurred earlier for alleles associated with low-level resistance. Understanding factors influencing acquisition and evolution of FQ resistance could reveal strategies for improved treatment success.

An Early Pandemic Analysis of SARS-CoV-2 Population Structure and Dynamics in Arizona.

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China, causing severe morbidity and mortality. Since then, the virus has swept across the globe, causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, the Arizona Department of Health Services, and those collected as part of community surveillance projects at Arizona State University and the University of Arizona. Phylogenetic analysis of 84 genomes from across Arizona revealed a minimum of 11 distinct introductions inferred to have occurred during February and March. We show that >80% of our sequences descend from strains that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related case in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020.IMPORTANCE As the COVID-19 pandemic swept across the United States, there was great differential impact on local and regional communities. One of the earliest and hardest hit regions was in New York, while at the same time Arizona (for example) had low incidence. That situation has changed dramatically, with Arizona now having the highest rate of disease increase in the country. Understanding the roots of the pandemic during the initial months is essential as the pandemic continues and reaches new heights. Genomic analysis and phylogenetic modeling of SARS-COV-2 in Arizona can help to reconstruct population composition and predict the earliest undetected introductions. This foundational work represents the basis for future analysis and understanding as the pandemic continues.

Genomic variant-identification methods may alter Mycobacterium tuberculosis transmission inferences.

Pathogen genomic data are increasingly used to characterize global and local transmission patterns of important human pathogens and to inform public health interventions. Yet, there is no current consensus on how to measure genomic variation. To test the effect of the variant-identification approach on transmission inferences for Mycobacterium tuberculosis, we conducted an experiment in which five genomic epidemiology groups applied variant-identification pipelines to the same outbreak sequence data. We compared the variants identified by each group in addition to transmission and phylogenetic inferences made with each variant set. To measure the performance of commonly used variant-identification tools, we simulated an outbreak. We compared the performance of three mapping algorithms, five variant callers and two variant filters in recovering true outbreak variants. Finally, we investigated the effect of applying increasingly stringent filters on transmission inferences and phylogenies. We found that variant-calling approaches used by different groups do not recover consistent sets of variants, which can lead to conflicting transmission inferences. Further, performance in recovering true variation varied widely across approaches. While no single variant-identification approach outperforms others in both recovering true genome-wide and outbreak-level variation, variant-identification algorithms calibrated upon real sequence data or that incorporate local reassembly outperform others in recovering true pairwise differences between isolates. The choice of variant filters contributed to extensive differences across pipelines, and applying increasingly stringent filters rapidly eroded the accuracy of transmission inferences and quality of phylogenies reconstructed from outbreak variation. Commonly used approaches to identify M. tuberculosis genomic variation have variable performance, particularly when predicting potential transmission links from pairwise genetic distances. Phylogenetic reconstruction may be improved by less stringent variant filtering. Approaches that improve variant identification in repetitive, hypervariable regions, such as long-read assemblies, may improve transmission inference.

Rhizopus microsporus Infections Associated with Surgical Procedures, Argentina, 2006-2014.

Rhizopus spp. fungi are ubiquitous in the environment and a rare but substantial cause of infection in immunosuppressed persons and surgery patients. During 2005-2017, an abnormally high number of Rhizopus infections in surgery patients, with no apparent epidemiologic links, were reported in Argentina. To determine the likelihood of a common source of the cluster, we performed whole-genome sequencing on samples collected during 2006-2014. Most isolates were separated by >60 single-nucleotide polymorphisms, and we found no evidence for recombination or nonneutral mutation accumulation; these findings do not support common source or patient-to-patient transmission. Assembled genomes of most isolates were ≈25 Mbp, and multiple isolates had substantially larger assembled genomes (43-51 Mbp), indicative of infections with strain types that underwent genome expansion. Whole-genome sequencing has become an essential tool for studying epidemiology of fungal infections. Less discriminatory techniques may miss true relationships, possibly resulting in inappropriate attribution of point source.

Isolation of Anti-Inflammatory and Epithelium Reinforcing Bacteroides and Parabacteroides Spp. from A Healthy Fecal Donor.

Altered intestinal microbiota is associated with systemic and intestinal diseases, such as inflammatory bowel disease (IBD). Dysbiotic microbiota with enhanced proinflammatory capacity is characterized by depletion of anaerobic commensals, increased proportion of facultatively anaerobic bacteria, as well as reduced diversity and stability. In this study, we developed a high-throughput in vitro screening assay to isolate intestinal commensal bacteria with anti-inflammatory capacity from a healthy fecal microbiota transplantation donor. Freshly isolated gut bacteria were screened for their capacity to attenuate Escherichia coli lipopolysaccharide (LPS)-induced interleukin 8 (IL-8) release from HT-29 cells. The screen yielded a number of Bacteroides and Parabacteroides isolates, which were identified as P. distasonis, B. caccae, B. intestinalis, B. uniformis, B. fragilis, B. vulgatus and B. ovatus using whole genome sequencing. We observed that a cell-cell contact with the epithelium was not necessary to alleviate in vitro inflammation as spent culture media from the isolates were also effective and the anti-inflammatory action did not correlate with the enterocyte adherence capacity of the isolates. The anti-inflammatory isolates also exerted enterocyte monolayer reinforcing action and lacked essential genes to synthetize hexa-acylated, proinflammatory lipid A, part of LPS. Yet, the anti-inflammatory effector molecules remain to be identified. The Bacteroides strains isolated and characterized in this study have potential to be used as so-called next-generation probiotics.

Genomic characterization of Parengyodontium americanum sp. nov.

Modern genome analysis and phylogenomic methods have increased the number of fungal species, as well as enhanced appreciation of the degree of diversity within the fungal kingdom. In this context, we describe a new Parengyodontium species, P. americanum, which is phylogenetically related to the opportunistic human fungal pathogen P. album. Five unusual fungal isolates were recovered from five unique and confirmed coccidioidomycosis patients, and these isolates were subsequently submitted to detailed molecular and morphological identification procedures to determine identity. Molecular and morphological diagnostic analyses showed that the isolates belong to the Cordycipitaceae. Subsequently, three representative genomes were sequenced and annotated, and a new species, P. americanum, was identified. Using various genomic analyses, gene family expansions related to novel compounds and potential for ability to grow in diverse habitats are predicted. A general description of the genomic composition of this newly described species and comparison of genome content with Beauveria bassiana, Isaria fumosorosea and Cordyceps militaris shows a shared core genome of 6371 genes, and 148 genes that appear to be specific for P. americanum. This work provides the framework for future investigations of this interesting fungal species.

Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity.

The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2's unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.

Whole-genome and targeted sequencing of drug-resistant Mycobacterium tuberculosis on the iSeq100 and MiSeq: A performance, ease-of-use, and cost evaluation.

BACKGROUND: Accurate, comprehensive, and timely detection of drug-resistant tuberculosis (TB) is essential to inform patient treatment and enable public health surveillance. This is crucial for effective control of TB globally. Whole-genome sequencing (WGS) and targeted next-generation sequencing (NGS) approaches have potential as rapid in vitro diagnostics (IVDs), but the complexity of workflows, interpretation of results, high costs, and vulnerability of instrumentation have been barriers to broad uptake outside of reference laboratories, especially in low- and middle-income countries. A new, solid-state, tabletop sequencing instrument, Illumina iSeq100, has the potential to decentralize NGS for individual patient care. METHODS AND FINDINGS: In this study, we evaluated WGS and targeted NGS for TB on both the new iSeq100 and the widely used MiSeq (both manufactured by Illumina) and compared sequencing performance, costs, and usability. We utilized DNA libraries produced from Mycobacterium tuberculosis clinical isolates for the evaluation. We conducted WGS on three strains and observed equivalent uniform genome coverage with both platforms and found the depth of coverage obtained was consistent with the expected data output. Utilizing the standardized, cloud-based ReSeqTB bioinformatics pipeline for variant analysis, we found the two platforms to have 94.0% (CI 93.1%-94.8%) agreement, in comparison to 97.6% (CI 97%-98.1%) agreement for the same libraries on two MiSeq instruments. For the targeted NGS approach, 46 M. tuberculosis-specific amplicon libraries had 99.6% (CI 98.0%-99.9%) agreement between the iSeq100 and MiSeq data sets in drug resistance-associated SNPs. The upfront capital costs are almost 5-fold lower for the iSeq100 ($19,900 USD) platform in comparison to the MiSeq ($99,000 USD); however, because of difference in the batching capabilities, the price per sample for WGS was higher on the iSeq100. For WGS of M. tuberculosis at the minimum depth of coverage of 30x, the cost per sample on the iSeq100 was $69.44 USD versus $28.21 USD on the MiSeq, assuming a 2 × 150 bp run on a v3 kit. In terms of ease of use, the sequencing workflow of iSeq100 has been optimized to only require 27 minutes total of hands-on time pre- and post-run, and the maintenance is simplified by a single-use cartridge-based fluidic system. As these are the first sequencing attempts on the iSeq100 for M. tuberculosis, the sequencing pool loading concentration still needs optimization, which will affect sequencing error and depth of coverage. Additionally, the costs are based on current equipment and reagent costs, which are subject to change. CONCLUSIONS: The iSeq100 instrument is capable of running existing TB WGS and targeted NGS library preparations with comparable accuracy to the MiSeq. The iSeq100 has reduced sequencing workflow hands-on time and is able to deliver sequencing results in <24 hours. Reduced capital and maintenance costs and lower-throughput capabilities also give the iSeq100 an advantage over MiSeq in settings of individualized care but not in high-throughput settings such as reference laboratories, where sample batching can be optimized to minimize cost at the expense of workflow complexity and time.

A global to local genomics analysis of Clostridioides difficile ST1/RT027 identifies cryptic transmission events in a northern Arizona healthcare network.

Clostridioides difficile is a ubiquitous, diarrhoeagenic pathogen often associated with healthcare-acquired infections that can cause a range of symptoms from mild, self-limiting disease to toxic megacolon and death. Since the early 2000s, a large proportion of C. difficile cases have been attributed to the ribotype 027 (RT027) lineage, which is associated with sequence type 1 (ST1) in the C. difficile multilocus sequence typing scheme. The spread of ST1 has been attributed, in part, to resistance to fluoroquinolones used to treat unrelated infections, which creates conditions ideal for C. difficile colonization and proliferation. In this study, we analysed 27 isolates from a healthcare network in northern Arizona, USA, and 1352 publicly available ST1 genomes to place locally sampled isolates into a global context. Whole genome, single nucleotide polymorphism analysis demonstrated that at least six separate introductions of ST1 were observed in healthcare facilities in northern Arizona over an 18-month sampling period. A reconstruction of transmission networks identified potential nosocomial transmission of isolates, which were only identified via whole genome sequence analysis. Antibiotic resistance heterogeneity was observed among ST1 genomes, including variability in resistance profiles among locally sampled ST1 isolates. To investigate why ST1 genomes are so common globally and in northern Arizona, we compared all high-quality C. difficile genomes and identified that ST1 genomes have gained and lost a number of genomic regions compared to all other C. difficile genomes; analyses of other toxigenic C. difficile sequence types demonstrate that this loss may be anomalous and could be related to niche specialization. These results suggest that a combination of antimicrobial resistance and gain and loss of specific genes may explain the prominent association of this sequence type with C. difficile infection cases worldwide. The degree of genetic variability in ST1 suggests that classifying all ST1 genomes into a quinolone-resistant hypervirulent clone category may not be appropriate. Whole genome sequencing of clinical C. difficile isolates provides a high-resolution surveillance strategy for monitoring persistence and transmission of C. difficile and for assessing the performance of infection prevention and control strategies.

Genomic Analyses of Acute Flaccid Myelitis Cases among a Cluster in Arizona Provide Further Evidence of Enterovirus D68 Role.

Enteroviruses are a common cause of respiratory and gastrointestinal illness, and multiple subtypes, including poliovirus, can cause neurologic disease. In recent years, enterovirus D68 (EV-D68) has been associated with serious neurologic illnesses, including acute flaccid myelitis (AFM), frequently preceded by respiratory disease. A cluster of 11 suspect cases of pediatric AFM was identified in September 2016 in Phoenix, AZ. To determine if these cases were associated with EV-D68, we performed multiple genomic analyses of nasopharyngeal (NP) swabs and cerebrospinal fluid (CSF) material from the patients, including real-time PCR and amplicon sequencing targeting the EV-D68 VP1 gene and unbiased microbiome and metagenomic sequencing. Four of the 11 patients were classified as confirmed cases of AFM, and an additional case was classified as probable AFM. Real-time PCR and amplicon sequencing detected EV-D68 virus RNA in the three AFM patients from which NP swabs were collected, as well as in a fourth patient diagnosed with acute disseminated encephalomyelitis, a disease that commonly follows bacterial or viral infections, including enterovirus. No other obvious etiological causes for AFM were identified by 16S or RNA and DNA metagenomic sequencing in these cases, strengthening the likelihood that EV-D68 is an etiological factor. Herpes simplex viral DNA was detected in the CSF of the fourth case of AFM and in one additional suspect case from the cluster. Multiple genomic techniques, such as those described here, can be used to diagnose patients with suspected EV-D68 respiratory illness, to aid in AFM diagnosis, and for future EV-D68 surveillance and epidemiology.IMPORTANCE Enteroviruses frequently result in respiratory and gastrointestinal illness; however, multiple subtypes, including poliovirus, can cause severe neurologic disease. Recent biennial increases (i.e., 2014, 2016, and 2018) in cases of non-polio acute flaccid paralysis have led to speculations that other enteroviruses, specifically enterovirus D68 (EV-D68), are emerging to fill the niche that was left from poliovirus eradication. A cluster of 11 suspect cases of pediatric acute flaccid myelitis (AFM) was identified in 2016 in Phoenix, AZ. Multiple genomic analyses identified the presence of EV-D68 in the majority of clinical AFM cases. Beyond limited detection of herpesvirus, no other likely etiologies were found in the cluster. These findings strengthen the likelihood that EV-D68 is a cause of AFM and show that the rapid molecular assays developed for this study are useful for investigations of AFM and EV-D68.